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ABSTRACT 

The aim of this work is to propose a joint exploitation of heterogeneous datasets 
from high-resolution/few-channel experiments and low-resolution/many-channel ex- 
periments by using a multiscale needlet Internal Linear Combination (ILC), in order 
to optimize the thermal Sunyaev-Zeldovich (SZ) effect reconstruction at high resolu- 
tion. We highlight that needlet ILC is a powerful and tunable component separation 
method which can easily deal with multiple experiments with various specifications. 
Such a multiscale analysis renders possible the joint exploitation of high-resolution 
and low-resolution data, by performing for each needlet scale a combination of some 
specific channels, either from one dataset or both datasets, selected for their relevance 
to the angular scale considered, thus allowing to simultaneously extract high resolu- 
tion SZ signal from compact clusters and remove Galactic foreground contamination 
at large scales. 

Key words: methods: data analysis - galaxies: clusters: general - cosmic background 
radiation 



1 INTRODUCTION 

Multifrequency observations of the microwave sky are a mix- 



to detect new galaxy clusters at any redshift (Birkinshaw 



ture of various diffuse and compact components (Bouchet 



et al. 



1995 



Bouchet &; Gispert 1999): the Cosmic Mi- 



crowave Background (CMB) emission, the SZ effect from 
galaxy clusters, the foreground emission from the Galactic 
interstellar medium (ISM), the Cosmic Infrared Background 
(CIB) emission, the emission from compact extra-galactic 
radio sources, the instrumental noise. Each of these com- 
ponents has a distinctive frequency signature (which can 
be known or not), a distinctive spatial distribution on the 
celestial sphere, and a distinctive spectral distribution on 
the angular scales (power spectrum). The separation of the 
sky components in such multifrequency observations of the 
CMB is an important part of the processing and analysis 
of such observational data. The thermal Sunyaev-Zeldovich 
(TSZ) effect from galaxy clusters is a spectral distortion 
of the CMB black body radiation due to inverse Compton 
scattering of the CMB photons off hot electrons contained 



in the intra-cluster gas (Sunyaev & Zeldovich 1972| . It is 
responsible for secondary temperature anisotropics, which 
introduce an excess of power in the primordial CMB tem- 
perature anisotropics at small angular scales. This spectral 
distortion is independent of the cosmological redshift so that 
the measure of the TSZ effect is a powerful and unique tool 
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1999 Carlstrom, Holder & Reese 2002). The separation and 



the extraction of the TSZ effect from the other sky emis- 
sions is made possible by a multifrequency coherent analy- 
sis because of the prior knowledge of the frequency signa- 
ture of the SZ effect. The main limitation in detecting SZ 
clusters comes from the achievable resolution of the instru- 
ment and the level of contamination by the other sky emis- 
sions and the instrumental noise. Various component sepa- 
ration methods have been developed to extract the thermal 
SZ signal from the observed frequency maps of a single ex- 
periment. Some of them require a prior assumption on the 
template SZ profile such as the matched filtering methods 
Haehnelt fe Tegmark| (|1996|), |Herranz et al.| ( |2002a|b| ), 



and Melin, Bartlett & Delabrouille (2006). Other methods 



are blind such as the methods based on statistical indepen- 
dence (ICA in|Hyvarinen](|1999|);|Maino et al.| (|2002|), Spec- 



tral Matching ICA in |Delabrouille, Cardoso Sz Patanchon| 
fl2003| ); [Cardoso et al.| ( 2008| ), GMCA in |Bobin et al.| | |2008f 
which also benefits from the sparsity of the components). 
Some methods are parametric and require physical mod- 



eling of some components (MEM in Hobson et al. (1998), 
Commander in |Eriksen et al.| ( |2008| )). Non-parametric meth- 
ods include needlet Internal Linear Combination (NILC) in 
Leach et al.l (12008} and inlDelabrouille et al.M2009L MILCA 



in Hurier, Hildebrandt & Macias- Perez (2010), multidimen- 



sional NILC in Remazeilles, Delabrouille & Cardoso (2011a) 



and Remazeilles, Delabrouille k, Cardoso (2011b)). The ILC, 
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which has been used first on the data of the WMAP mission 

is 



2003) 



to reconstruct the CMB emission (Bennett et al. 
a particularly simple component separation method which 
does not assume any particular parametrization for fore- 
ground emission, and for the other emissions considered as 
contaminants. The frequency scaling vector of the compo- 
nent of interest to extract is the only parameter needed for 
implementing ILC. For all the component separation meth- 
ods, the problem is to find the best compromise between 
the most accurate unbiased estimation of the component of 
interest and the best minimization of the residuals from the 
other sky emissions. 



In this work we address the problem of jointly exploit- 
ing the datasets of multiple experiments with various spec- 
ifications for optimizing the reconstruction of the thermal 
SZ effect at high resolution. No single instrument provides 
the best measurement of the SZ effect. Space-borne instru- 



ments, such as Planck for instance ( 


Planck Collaboration 


2011a||PPlanck Collaboration et al.| 


2011| |Planck Collabo-| 


ration et al.| 2012b}, can observe the whole sky in the mil- 



limetre through many frequencies (nine channels from 30 
GHz to 857 GHz). The large number of frequency channels 
covered by such space-borne instruments enables component 
separation methods, such as ILC, to improve the minimiza- 
tion of various contaminations (Galactic foregrounds, CMB, 
instrumental noise) in the reconstructed SZ map. However, 
the instrumental resolution of Planck for instance is limited 
to 5 arcmin, therefore the extraction of SZ signal from clus- 
ters of arcminute angular size, by using such dataset only, 
remains unachievable. Conversely, higher resolution ground- 
based telescopes, such as the Atacama Cosmology Telescope 
(ACT) (|Marriage et al.||2011|) and the South Pole Telescope 



(SPT) (Williamson et al. 2011 Reichardt et al. 2012) for 



instance, are designed to measure the SZ temperature an- 
isotropics up to arcminute angular scales but the presence 
of the atmosphere limits the frequency coverage of those ex- 
periments. The few number of frequency channels observed 
by such ground-based telescopes (148 GHz, 218 GHz, and 
277 GHz in ACT) limits the ability of these instruments 
to remove the Galactic foreground contamination in the re- 
construction of the SZ signal of nearby clusters. Here we 
propose to perform SZ component separation from the joint 
exploitation of heterogeneous sets of maps (both maps from 
many-channel/lower-resolution instrument and maps from 
few-channel/higher-resolution instrument) by using needlet 
ILC as a multiscale, or multiresolution, approach for com- 
bining multiple instrument datasets. 



The paper is organised as follows. In Sect. [2] we de- 
scribe the needlet ILC method on patches of the sky and 
show how the needlets can be exploited to combine multiple 
instrument datasets with various resolutions for optimizing 
the SZ component separation. The results are presented on 
simulations of the sky in Sect. [3] where we apply the method 
jointly on two heterogeneous instrument datasets. We dis- 
cuss the robustness of the approach in Sect. ^] and conclude 
m Sect. [5] 



2 AGGREGATING DATASETS WITH 
NEEDLET ILC 

2.1 Derivation of ILC 

All the available observation maps (N Q h s maps, indexed by 
i for each frequency channel) can be written as 



Xi{p) = ais(p) +rii(p), 

which can be recast in a vector form as 

x(p) = as(p) +n(p). 



(i) 



(2) 



Here p may denote either a pixel in a direct-space represen- 
tation, or an (£,m) index (resp. k — (ki, fe)) in a harmonic 
space (resp. Fourier space) representation, or even a wavelet 
domain (j;q) in a wavelet frame decomposition, where j 
denotes the scale of the wavelet and q the pixel. The vec- 
tor x(p) collects the iV bs observed frequency maps, s(p) is 
the unknown thermal SZ template map that we would like 
to reconstruct, and a is the known frequency scaling vec- 
tor of the SZ effect. All the other sky emissions (CMB and 
foregrounds) and the instrumental noise are collected in a 
single nuisance term n(p). The ILC estimate of the SZ map, 
*s — WiXi(p), is a linear combination of the observed maps 
of minimum variance under the constraint of offering unit 
response to the SZ component, so that the ILC weights w 
are solution to the following constrained minimum variance 
problem: 



(3) 



J var(?(p)) minimum, 

\ WiCLi = 1 

The first condition in Eq. ^ guarantees the minimum con- 
tamination by the background noise (residual sky compo- 
nents and instrumental noise) whereas the second condition 
guarantees the unbiased reconstruction of the SZ template. 
Note that the quality of the reconstruction relies on the accu- 
rate knowledge of the SZ frequency scaling a. In presence of 
calibration errors, there is no guarantee that the SZ compo- 
nent is conserved (Dick, Remazeilles & Delabrouille 2010). 



Straightforward algebra using Lagrange multiplier (Eriksen 
et al.| 2004| ) leads to the following ILC estimate as the solu- 
tion to Eq. 



s(p) 



r x(p) 



a 1 R 



-x(p), 



(4) 



where R^/ = -^j-J2 P x i(p) x i'(p) * s the empirical covari- 
ance matrix of the observations. By construction, we have 
s"(p) = s(p) + w T n(p), with the residual nuisance term min- 
imized. The ILC relies on the component of interest to be 
uncorrelated with the contaminants (i.e. (s(p)rii(p)) = 0, 
for each channel). We discuss the case where this assump- 
tion can fail in Sect. |4.l] when applying ILC on small patches 
of the sky. 

Let us now discuss the drawbacks of using either pixel- 
based or multipole-based ILC filtering in the context of ex- 
ploiting heterogeneous datasets from multiple instruments 
with different specifications. 



The ILC filtering in pixel space (Eriksen et al. 2004), 
which is ultra-local in direct space, suffers from nonlocal- 
ity in harmonic (or Fourier) space by conjugation. Such a 
pixel-based ILC thus requires to combine the same set of 
channel-maps for all the angular scales considered. This is 
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not optimal when the channel-maps have not the same an- 
gular resolution since it requires that the fixed set of maps 
be degraded to a common resolution which is compelled to 
be the lowest resolution of the set of channel-maps, so that 
the ILC reconstructed map is also compelled to be produced 
at this lowest resolution and the small scale information is 
lost. 

Consequently, pixel-based ILC is either a resolution- 
limited method or a channel-limited method: suppose that 
we would like to reconstruct an SZ map at high resolution, 
say one arcminute, we would have to combine only the few 
channel-maps having a resolution which is higher than or 
equal to one arcminute. A reduced number of channels for 
SZ ILC reconstruction is not optimal for removing the con- 
tamination by other sky components because the number of 
sky components becomes larger than the number of chan- 
nels, making the inverse problem suboptimal (also see the 
results in Sect. [3|. Conversely, if we choose to combine a 
larger number of channel-maps, for instance by including in 
the pixel-based ILC combination lower resolution channel- 
maps from another instrument, in order to better remove 
the background noise, then the ILC output SZ map would 
be restricted to have the lowest original resolution of the 
combined channel-maps. 

As an alternative, the component separation can in 
principle be performed in harmo nic space (|Tegmark, de 
' Oliveira-Costa fe Hamilton) [20031 l Kim ? Naselsky fe Chrls^ 
tensen[ |2008| ). A multipole-based ILC can mix maps with 
different resolutions, unlike a pixel-based method, by com- 
bining all the channel-maps at low £ and only the high res- 
olution channel-maps at higher L A harmonic space ILC 
(ultra-local in £) is a fortiori nonlocal in pixel space so that 
such an approach is better suited for all-sky analysis. There- 
fore, a multipole-based ILC can not easily handle data from 
a partial sky coverage. Ground-based experiments cover fi- 
nite size sky area. The advantage of such partial sky surveys 
is that one is left with analyzing zones of the sky which are 
less contaminated by other sky emissions such as Galactic 
foregrounds. Here we are interested on performing compo- 
nent separation on patches of the sky thus we cannot use a 
harmonic space approach. 

2.2 Needlet ILC: mult iresolut ion 

For the purposes of a joint analysis of heterogeneous datasets 
(many channels- lower resolution/few channels-higher resolu- 
tion) , we propose to use a wavelet-based component separa- 



loux, Fay & Cardoso ( 2009 ) . The needlet decomposition pro- 



tion called Needlet ILC ( Delabrouille et al. 2009) because 



of the ability of a wavelet-based approach to combine data 
with different resolutions, unlike a pixel-based approach, and 
to handle patches of the sky, unlike a harmonic space ap- 
proach. The needlets are a particular type of wavelets where 
the widths of the scale windows are freely tunable in har- 



monic space (Guilloux, Fay & Cardoso 2009). The ability of 



wavelets for combining multi-instrument datasets has been 
highlighted in the case of CMB power spectrum reconstruc- 
tion in |Fay et aL] ( |2008 ). 

Instead of requiring strict locality either in pixel or in 
multipole space, we can perform a needlet-based component 
separation by decomposing the maps on a needlet frame. 
The formulation of the needlets in the context of CMB data 



analysis can be found in Marinucci et al. (2008) and GuiL 



vides localization of the ILC filter both in pixel and in har- 
monic space, making the ILC filtering able to adapt to the 
local conditions of contamination (local Galactic contami- 
nation, small scale noise contamination). Wavelet localiza- 
tion has also been adopted by Bobin et al. (2012) in the 



GMCA method, based on sparsity, for optimizing compo- 
nent separation. The property of localization of the needlets 
both in scale and in space is particularly meaningful for 
an SZ dedicated component separation because the signal 
is non-uniformly distributed both in space (on the celestial 
sphere) and in scale (relevant at small scales). Furthermore, 
Needlet ILC (NILC) is a multiresolution approach allow- 
ing for combining channel-maps with different resolutions 
by using localization in Fourier or harmonic space (needlet 
bands). Moreover, such a wavelet approach allows us to com- 
bine patches of the sky by using localization in pixel space. 
Let us define a collection of window functions (Fig. [T]) 
in the Fourier domain, indexed by j, such that over the 
useful range of spatial frequency k, we have 



(5) 



In the following we will perform component separation on 
patches of the sky smaller than 100 square degrees so that 
we adopt a flat-sky approximation allowing for using Fast 
Fourier transforms (FFT) instead of spherical harmonic 
transforms. Maps of needlet coefficients x^ j \p) at a given 
scale (j) are obtained, for the observed maps x(p), by inverse 
FFT of the associated map Fourier coefficients Xk filtered 
by the spectral windows hf^' 



N k -lN k - 



k 1 =0 k 2 =0 



i2n 



ki , k 2 



(6) 



where p = (pi, ^2), k = \Jk\ + k\ and Nk is the number of 
resolution elements on each side of the two-dimensional FFT 
patch considered. The Needlet ILC estimate at the scale (j) 
is given by 



s U) (p) = w T x U) (p) 



aT(R U \p))-ia 



x U) (p), 



where the local covariance matrix of the needlet coefficients 
of the observations is 



(8) 



The pixel domain T> p , on which the local covariance is com- 
puted, is defined from the smoothing of the product map 
xi j \p)x[ j \p) with a symmetric Gaussian window in pixel 
space. The complete reconstruction of the SZ map s~(p) can 
then be obtained by simply coadding the reconstructed SZ 
"maps per scale" z^ j \p), 



(9) 



where the SZ "maps per scale" themselves are obtained from 
the following scheme 



Sk 



\p)- 



(10) 
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Figure 1. Illustration of needlet scales versus beams. In this plot 
we assume N\ = 5 and N2 = 4. The beam of expi channels 
(thick dashed red) and the beam of EXP2 channels (thick solid 
blue) are overplotted on the needlet bands (thin dotted black). For 
each needlet band an ILC estimate is computed: in the first four 
bands the NILC processing can combine expi and EXP2 channels 
whereas in the last four bands only EXP2 channels are exploited 
by NILC. 



The constraint on the needlet scale-bands in Eq. |5| guaran- 
tees the total conservation of the SZ power by the processing 
in Eqs. $\ and (To). 



2.3 



Example of needlets for heterogeneous 
datasets 



The needlet ILC component separation has already been 
dedicated to all-sky TSZ reconstruction ( |Leach et al .[|2008). 



Here we develop needlet ILC on patches of sky to focus on 
zones of the sky centered on particular clusters. Our main 
purpose is to use needlets on patches as a multiresolution 
tool for combining multiple heterogeneous experiments by 
ILC to optimize the TSZ component separation. 

Let us illustrate how such a multiscale approach to com- 
ponent separation is relevant to the joint exploitation of het- 
erogeneous data with different resolutions. In this regard we 
consider two experiments noted expi and EXP2, the former 
having a lower resolution 0\ than the latter having a beam 
O2 < Oi , but a larger number of frequency channels Ni > N2 . 
For the sake of simplicity we assume that the Ni (resp. iV 2 ) 
channel-maps of expi (resp. EXP2) have the same resolution 
61 (resp. 62). The goal is to benefit both from the maximum 
number of channels Ni + N2 and from the highest resolu- 
tion 62 for performing NILC component separation, both re- 
quirements being needed to guarantee a reconstructed TSZ 
signal with the lowest sky residual contamination and with 
the highest resolution. 

We may define Ni + N2 needlet bands such that 



the beam of expi, #1, can not probe the scales k > k : 



where k^L^ is the highest Fourier mode probed by the A^ h 
needlet band h[ Nl) (Fig.[l]). 

It is then possible to perform ILC filtering at each 
needlet scale (j) by combining the channel-maps whose res- 
olutions are compatible with the needlet band consid- 
ered. Typically, in the bands (j), where 1 < j < N±, the ILC 



TSZ output (p) is obtained from the combination of the 
N1+N2 maps of both expi and EXP2. Next, in the following 
bands (j), where Ni + 1 ^ j ^ N2, the ILC combination is 
performed on the N2 higher resolution maps of EXP2 only, 
the maps of expi having no information at the scales cov- 
ered by the N2 latest bands to be taken into account by the 
ILC combination. Formally, for all j ^ Ni the ILC estimate 
is computed from the combination of both the expi dataset 
and the EXP2 dataset 



(p) 



U) Aj) 



+ ^nWnI bxpx(p) 

n 0') 



U N 1 +1^N 1 + 1 exp 2 (p) + ••• + Wn[+n 2 X N\+N 2 exp 

while for all j > Ni the ILC is computed from the combi- 
nation of the EXP2 dataset only 



SO") 



(p) 



„(i) 



S3) 



» + - + 



iVi + l^iVi + 1 EXP2V-P; ~ ••• ~ Ni IJY2 JVj +A/"2 exp 2 

Therefore, NILC appears as a tunable multiresolution 
tool to aggregate multiple experiments with various specifi- 
cations: larger angular scales are exploited to remove con- 
taminants (Galactic foregrounds, etc) from the NILC combi- 
nation of a large number of channels (both expi and EXP2), 
whereas smaller angular scales are exploited to reconstruct 
SZ effect at higher resolution from the combination of the 
high resolution EXP2 channels only. 



3 RESULTS 

3.1 Sky simulations 

We now illustrate our method on the high r esolution sim- 
ulations of the microwave sk}|^] generated by |Sehgal et al. 
(2010). For this study the simulated maps have been de- 
graded from their original HEALPix resolution nside=8192 
to a HEALPix resolution nside=4096. The simulated all-sky 
emission is generated along six frequency channels (30 GHz, 
90 GHz, 148 GHz, 219 GHz, 277 GHz and 350 GHz). We 
have made use of the CMB emission, the thermal SZ effect, 
the Galactic dust emission, and the emissions from infrared 
and radio galaxies. All these components are coadded in each 
frequency channel to provide six sky maps. From these sim- 
ulations we have created two datasets, a Planck-like set of 
maps and an ACT-like set of maps. The Planck-like and the 
ACT-like experimental characteristics are defined in Tab.[l] 
Planck-like maps are constructed by smoothing the simu- 
lated sky maps, except the 277 GHz sky map, to the beams 
given in Tab. [l] and by adding a white noise map from the 
RMS levels of Tab. [l] At this stage we have five Planck- like 
channel-maps but we also created a sixth map by extrap- 
olating the 350 GHz simulated components to 545 GHz. 
We have used for this extrapolation the temperature and 
spectral index of Galactic dust computed in Planck Collab- 
(2011c|, and the temperature and spectral index of 



oration 



infrared sources computed in Planck Collaboration (2011b) 



ACT-like maps are constructed by smoothing the 148 GHz, 
219 GHz, and 277 GHz simulated sky maps to the beams 
given in Tab. j^and by adding a white noise map at RMS 
values of Tab.^l 

In the next sections, we apply NILC on three sets 



The 



simulation 



be 



found 



at 



http:/ / lambda, gsfc. nasa.gov/toolbox/tb_cmbsim_ov.cfm 
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Planck-like 


30 GHz 


90 GHz 


148 GHz 


219 GHz 


350 GHz 


545 GHz 


[arcmin] 


32.65 


9.88 


7.18 


4.71 


4.5 


4.72 


<Tjv(4 surveys) 


210 


57 


30 


44 


173 


2485 


[fiK • arcmin] 















ACT-like 


148 GHz 


219 GHz 


277 GHz 


[arcmin] 


1.4 


1.0 


0.9 


<7jv [/^K • arcmin] 


30 


45 


60 



Table 1. Specifications of both datasets constructed from simulations o f|Sehgal et al.| ( |2010| ). Beams and noise RMS for Planck-like maps 
refer to the values in Mennell a et al.| |2011| ), [pfanck HFI Core Team| ( |2011| > extrapolated to 4 full-sky surveys (28 months). ACT-like 
values refer to the 2008 ACT Southern survey pas et al.||2011^ 




Figure 2. A 6.8° x 6.8° patch of the simulated sky centered on the selected cluster 1 of Tab. [2] On the top panels: input thermal SZ 
smoothed to 3 arcmin resolution (left) and to 10 arcmin resolution (right). On the bottom panels: simulated Galactic dust at 277GHz 
(left) and ACT-like 277GHz channel-map (right). 





z 


R.A. 


Dec. 


M 20 o 


^200 


#200 


^200 






[degree] 


[degree] 


[1O 14 M ] 


[Mpc] 


[arcmin] 


[arcmin 2 ] 


cluster 1 


0.04 


39.6 


15.8 


7.3 


1.82 


40.7 


0.0299 


cluster 2 


0.95 


2.5 


78.0 


9.7 


1.45 


3.0 


0.0007 



Table 2. Two selected clusters from the SZ Halo catalog of the simulation of Sehgal et al. (2010). For each cluster, the table collects its 
redshift, its right ascension (R.A.) and its declination (Dec), its mass, its radius, its angular size on the sky, and its integrated Compton 
value within R200. 




Figure 3. NILC thermal SZ reconstruction of the selected cluster 1 of Tab. [2] from different datasets at 10 arcmin (top panels) and 3 
arcmin (bottom panels). Left panels: NILC TSZ reconstruction from Planck-like channels only. Middle panels: NILC TSZ reconstruction 
from ACT-like channels only. Right panels: NILC TSZ reconstruction from the combination of heterogeneous data, i.e. both Planck-like 
and ACT-like channels. 



of data (Planck-like dataset, ACT-like dataset, and joint 
Planck- like/ ACT- like dataset) to retrieve the TSZ signal 
on 7° x 7° patches of the sky centered on a chosen clus- 
ter. The TSZ reconstruction will be performed at 10 and 5 
arcmin resolution^] when NILC is applied on the Planck- 
like dataset, whereas the reconstruction will be performed 
at 10 and 3 arcmin resolutions when NILC is applied either 
on the ACT-like dataset or on the joint Planck- like /ACT- 
like dataset. The reconstruction will be illustrated on two 
patches of the sky centered on two selected clusters given in 
Tab. [2] (an extended one with #200 = 40.7 arcmin and a com- 
pact one with #200 = 3 arcmin) and then generalised to a 
sample of selected clusters. Here we define #200 = R200/DA 
where R200 is the radius in which the total density contrast 
of the cluster is 5 = 200, as compared to the critical density 
of the Universe at the cluster redshift, and Da is the angular 
diameter distance to the cluster. 

3.2 NILC Thermal SZ reconstruction from 
combined ACT-like/Planck-like data 

The NILC reconstruction of the thermal SZ effect due to the 
selected extended cluster 1 of Tab. [2] with typical angular 
size #200 =40.7 arcmin and located at a redshift z = 0.04, is 
shown in the panels of Fig. [3] whereas the input signal from 
the simulations is shown in Fig. [2] We can see in the middle 

2 The smallest resolution that we can reach when NILC is applied 
on the single Planck-like dataset is about 5 arcmin (see Tab. (I). 



panels of Fig.[3]that the reconstructed TSZ map from ACT- 
like maps only is clearly more contaminated by foreground 
residuals (visible on the top right corner of the patches) than 
the one reconstructed from Planck-like maps only (left pan- 
els) . Combining both datasets for the reconstruction further 
improves the removal of foreground residuals (right panels 
of Fig. [3|, by taking advantage of the large number of chan- 
nels. Simultaneously, the combination of both Planck- like 
and ACT-like maps through Needlet ILC enables to recon- 
struct the thermal SZ signal at high resolution (3 arcmin 
in bottom right panel of Fig.j3|, which would be impossible 
from the use of Planck-like maps only since the beams are 
limited to 5 arcmin. Note that on the bottom left panel the 
reconstructed TSZ map from Planck-like data actually has 
a resolution of 5 arcmin. Moreover, Planck-like TSZ recon- 
struction at high resolution is noisy compared to ACT-like 
TSZ reconstruction because of the lower sensitivity. Once 
again, the joint exploitation of both datasets by NILC signif- 
icantly reduces both the instrumental noise contamination 
and the sky residuals (right panels of Fig. |3j). These results 
clearly show that a reconstruction of thermal SZ from the 
combination of heterogeneous datasets (with different num- 
ber of channels and different resolutions) benefits from the 
advantages of the TSZ reconstructions from each dataset in- 
dependently while avoiding their drawbacks: low foreground 
contamination from Planck-like channels, high resolution 
and low noise contamination from ACT-like channels. For 
illustration, the NILC "maps per scale", which provide the 
reconstructed TSZ map by coaddition, can be seen in Fig. 
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Figure 4. NILC reconstruction per scale (Planck-like/ ACT-like combined) for the selected cluster 1 of Tab. [2] Last panel: the associated 
needlet bands used for the reconstruction. 



[4] Each patch corresponds to the ILC reconstruction at a 
given range of angular scales defined by the needlet spectral 
bands plotted in the last panel of Fig. [4] 

The various power spectra of the reconstructed TSZ 
maps are shown in top panels of Fig. [5] The power spec- 
trum of the NILC reconstructed TSZ map from the joint 
exploitation of both Planck-like maps and ACT-like maps 
best reconstructs the input TSZ power spectrum for the 10 
arcmin reconstruction (top left panel of Fig. [5} whereas the 
3 arcmin reconstruction becomes noisy at small scales k > 
0.05 arcmin -1 . As expected, on the top right panel of Fig. [5] 
NILC on Planck-like data better reconstructs the TSZ power 
spectrum at large scales (k < 0.05 arcmin -1 ) than ACT- 
like NILC because of the larger number of available chan- 
nels used to remove the foreground contamination whereas 
NILC on ACT-like data performs better reconstruction at 
small scales than Planck-like NILC which lacks resolution 
and is noise-dominated. At large scales (k < 0.05 arcmin -1 ) 



NILC applied on joint datasets matches the Planck- like re- 
construction whereas at small scales (k > 0.05 arcmin -1 ) it 
matches the ACT-like reconstruction, thus performing the 
best compromise. By applying the NILC weights onto the 
simulated frequency dust maps we are able to estimate the 
contamination level of Galactic dust emission in the recon- 
structed TSZ maps from each dataset. The power spectra of 
the residual Galactic dust emission for the 3 arcmin recon- 
struction are plotted in the bottom left panel of Fig. [5] at 
scales k < 0.1 arcmin -1 , the ACT-like TSZ reconstruction 
shows more residual Galactic dust than the Planck-like TSZ 
reconstruction, as it is expected from the reduced number 
of exploited channels. At smaller scales k > 0.1 arcmin -1 , 
reconstruction in Planck-like case shows more Galactic dust 
residuals because the ILC weights are requested for minimiz- 
ing the instrumental noise which becomes the dominant con- 
taminant here, as it can be also seen in bottom right panel 
of Fig. [5] The TSZ NILC reconstruction from the joint ex- 
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Figure 5. NILC thermal SZ power spectra from different datasets for 10 arcmin reconstruction (top left panel) and 3 arcmin reconstruc- 
tion (top right panel) versus input power spectrum (selected cluster 1 of Tab. [2j. At the smallest scales, the Planck-like TSZ spectrum 
is noisy compared to the ACT-like TSZ power spectrum and limited to the instrumental beam. The joint Planck-like/ ACT-like TSZ 
spectrum matches very well the input thermal SZ spectrum at 10 arcmin resolution. Residual Galactic dust power spectra (bottom left 
panel) and residual instrumental noise power spectra (bottom right panel) for the 3 arcmin TSZ reconstruction are also plotted for 
comparison. At the largest scales, the spectrum of the reconstructed TSZ from Planck-like data is less contaminated by the residual 
Galactic dust than the TSZ power spectrum from ACT-like data because of the higher number of channels combined by ILC. The joint 
Planck-like/ ACT-like reconstruction shows an even better reduction of the residual dust power. 



ploitation of both ACT-like and Planck- like datasets reduces 
further the power of residual dust on all scales, matching the 
Planck-like reconstruction on large scales and the ACT-like 
reconstruction on small scales. In that case, on the largest 
scales, Galactic dust residuals are two orders of magnitude 
below the TSZ signal (bottom left panel of Fig. [5). The level 
of residual instrumental noise (bottom right panel of Fig. [5]) 
is also reduced in all scales when reconstructing on joint 
ACT-like and Planck-like datasets, unlike single instrument 
reconstructions. However, instrumental noise still dominates 
the TSZ power spectrum on scales k > 0.1 arcmin -1 . 

In Fig. [5] we plot the reconstructed TSZ cluster pro- 
files versus the input TSZ cluster profile for the NILC re- 
construction at 10 arcmin. The analogue, at 3 arcmin, is 
plotted in Fig. [7] We can see on the left panel the TSZ pro- 
file obtained from a reconstruction based on the Planck-like 
observations only, the middle panel shows the TSZ profile 
from a reconstruction based on the ACT-like observations 
only, and the right panel exhibits the profile from a recon- 
struction jointly exploiting both ACT-like and Planck-like 
channels. We see how the combination of complementary 
datasets in the NILC filtering improves the TSZ profile re- 



construction: less contamination is observed (right panels of 
Fig. [6]) compared to single instrument reconstructions (left 
and middle panels). Moreover, the NILC filtering on joint 
datasets is able to recover the TSZ profile at high resolution 
(right panels of Fig.[7| where NILC on Planck-like data only 
fails (left panels of Fig. We quantitatively measure the 
performance of the reconstruction of the SZ profile by the 
trade-off between the bias and the variance as a criterion. 
This trade-off is measured by the mean squared error (MSE) 

MSE = ^2(?i- y/ in) ) 2 + al (11) 

ring i 

where Y{ is the reconstructed profile, Y^ the input pro- 
file, and Gi the standard deviation (error bar) on the recon- 
structed profile for a ring i of pixels. The standard deviation 
is given by 

V P £ ring 

where y is the reconstructed SZ map and N p the number 
of pixels in the ring i considered. The results are listed in 
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Figure 6. NILC TSZ profiles of the selected cluster 1 of Tab. [2] (10 arcmin TSZ reconstruction). Planck-like (left panels), ACT-like 
(middle panels), Planck- like/ ACT- like combined (right panels). 
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Figure 7. NILC TSZ profiles of the selected cluster 1 of Tab. [2] (3 arcmin TSZ reconstruction). Planck-like (left panels), ACT-like 
(middle panels), Planck- like/ ACT- like combined (right panels). 



Tab. [3] for each reconstruction. We see from this criterion 
that Planck- like 10 arcmin reconstruction of the TSZ pro- 
file is better than ACT-like reconstruction whereas ACT-like 
NILC better reconstructs high resolution profile at 3 arcmin 
than Planck-like NILC. For both output resolutions, the best 
reconstruction of the TSZ profile is performed by the joint 
Planck-like/ ACT- like NILC. Table [3] indicates that Planck- 
like channels are exploited by NILC to remove sky contam- 
ination at large scales and lower resolution whereas high 
resolution ACT-like channels are exploited for reconstruct- 



ing the SZ signal at higher resolution. The oscillations in the 
reconstructed profile (Figs. [6] and [7| far from the center of 
the cluster (at a radius larger than 50 arcmin) are due to 
residual temperature fluctuations of sky contaminants and 
instrumental noise residuals. We can see on the bottom pan- 
els of Fig. [7] that the amplitude of oscillations is smaller 
when using both datasets for the reconstruction. Moreover, 
the typical wavelength of the oscillations is smaller in the 
Planck-like case than in the ACT-like case, indicating that 
Planck-like reconstruction is mostly contaminated by instru- 
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Figure 8. Residual dust profile (top panels) and residual noise profile (bottom panels) for the selected cluster 1 of Tab. [2] Left panels: 
10 arcmin TSZ reconstruction. Right panels: 3 arcmin TSZ reconstruction. Planck-like (dot-dashed black), ACT-like (long dashed red), 
Planck-like/ ACT-like combined (solid blue). 



MSE (10 arcmin) 


Planck- like 


ACT-like 


Joint 




1.77e-9 


3.82e-9 


1.20e-9 


MSE (3 arcmin) 


Planck- like 


ACT-like 


Joint 




2.70e-8 


9.58e-9 


5.11e-9 



Table 3. Performance of the profile reconstruction measured by 
the criterion Eq. ( |ll| >. Whereas Planck-like reconstruction is bet- 
ter than ACT-like reconstruction at 10 arcmin, it is the opposite 
for the 3 arcmin reconstruction. For both output resolutions, the 
best reconstruction is performed from the joint Planck-like/ ACT- 
like reconstruction. 



mental noise (small wavelength fluctuations) whereas ACT- 
like reconstruction is mostly contaminated by Galactic resid- 
uals (longer wavelength fluctuations). This is also confirmed 
by the right panels of Fig. [8] showing both residual Galac- 
tic dust in the profile reconstruction and residual noise with 
smaller wavelength fluctuations. These results illustrate how 
NILC applied jointly on Planck-like/ ACT-like datasets is 
able to reconstruct the TSZ effect with a low sky residual 
contamination on the one hand (Fig. |6| and at a resolution 
beyond the beam of the Planck-like instrument on the other 
hand (Fig. [7]). 



We also compute the contamination of Galactic dust 
and instrumental noise in the reconstructed TSZ profile by 
applying the NILC weights either onto the input dust maps 
of the simulation or onto the simulated noise maps. Top pan- 
els of Fig. [8] represent the contamination of Galactic dust 
in the reconstructed TSZ profile for the three sets of data 
(Planck- like only, ACT-like only, and joint datasets). Again 
we see that Galactic dust contaminates more the ACT-like 
TSZ profile reconstruction with larger amplitude fluctua- 
tions than the Planck-like one because of the reduced num- 
ber of exploited channels in the ACT-like case. The contam- 
ination of Galactic dust is clearly reduced further when per- 
forming NILC filtering jointly on both datasets. The residual 
instrumental noise in the TSZ profile on the bottom right 
panel of Fig. [8] shows oscillations with larger amplitude in 
the Planck-like case than in the ACT-like case for the 3 ar- 
cmin TSZ reconstruction, and the amplitude decreases when 
applying NILC jointly on both datasets. 

3.3 Thermal SZ reconstruction of a compact 
cluster 

We now repeat the analysis on a more compact cluster. The 
compact cluster 2 of Tab. [2] is located at a redshift z ~ 1 
with a typical angular size on the sky #200 — 3 arcmin. 

The reconstruction of the selected compact cluster 
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Figure 9. A 6.8° x 6.8° patch of the simulated sky centered on the selected compact cluster 2 of Tab. [2] A yellow circle has been 
overplotted to indicate the location of cluster. Top left panel: Input thermal SZ of the compact cluster smoothed to 3 arcmin resolution. 
Top right panel: NILC TSZ reconstruction from Planck-like channels only (limited to 5 arcmin reconstruction), no TSZ detection, 
instrumental noise contaminated. Bottom left panel: NILC TSZ reconstruction from ACT-like channels only (3 arcmin reconstruction), 
TSZ detection, sky contaminated. Bottom right panel: NILC TSZ reconstruction from the combination of both Planck-like and ACT-like 
channels (3 arcmin reconstruction), TSZ detection with low contamination. 



is shown in Fig. [9] Whereas an ILC applied on Planck- 
like channels only is unable to reconstruct the TSZ signal 
from such a cluster (top right panel of Fig. [9|, the same 
component separation applied on high-resolution ACT-like 
channel-maps successfully reconstructs the TSZ signal (bot- 
tom left panel of Fig. [9). However, the background residual 
noise from Galactic foregrounds is not successfully filtered. 
On the bottom right panel of Fig. [9] the wavelet ILC ex- 
ploiting both Planck-like and ACT-like channel-maps is able 
to simultaneously reconstruct the TSZ signal from such a 
compact cluster and to successfully clean the background 
residual noise. The joint exploitation of both datasets by 
a needlet ILC method shares with the ACT-like NILC the 
ability of reconstructing the TSZ effect from compact clus- 
ters of typical angular size below the beam of the single 
Planck-like instrument on the one hand, and shares with 
Planck-like NILC the ability of minimizing the sky residual 
contamination on the other hand. 



3.4 Statistical analysis on a sample of SZ galaxy 
clusters 

We now perform the same analysis on a selected sample of 90 



SZ clusters from the catalog of Sehgal et al. ( 2010 ) consisting 



in three sub-samples of 30 clusters each: a sub-sample of 30 
extended clusters with 15 arcmin < #200 < 40 arcmin, a sub- 
sample of 30 "typical" clusters with #200 00 5 arcmin, and 



a sub-sample of 30 compact clusters with #200 ~ 3 arcmin. 
Their different locations on the sky allow us to explore dif- 
ferent levels of sky background. The mean TSZ profiles are 
shown in Fig. |10| 

On the sub-sample of extended clusters (top panels of 



Fig. 10 ) the reconstruction of the TSZ profile by NILC is re- 



liable for the three sets of data (Planck-like only, ACT-like 
only, and joint datasets) with the ACT-like reconstruction 
slightly more contaminated by sky residual fluctuations at 
large radii. This trend is confirmed when we look at the con- 
tamination of Galactic dust in the reconstructed TSZ mean 
profile (Fig. 12): residual Galactic dust has a larger ampli- 
tude in the ACT-like reconstruction. When we are interested 
in reconstructing the TSZ profiles of compact clusters with 
#200 < 5 arcmin (middle panels and bottom panels of Fig. 
10) we can take advantage of the high resolution ACT-like 
channel-maps in the NILC combination to go beyond the 
limited beam of the Planck- like instrument. The middle left 
and bottom left panels of Fig. |10| confirm that a component 
separation applied on the single Planck-like dataset fails to 
recover the TSZ effect from compact clusters. 

For completeness we compute the mean squared error 
measuring the trade-off between the bias and the variance of 
the profile reconstruction by generalising the criterion Eq. 
( 11 ) for a single cluster to the case of a cluster sample: 
2 



MSE 



(13) 
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Figure 10. Cluster mean y-profiles for 7° x 7° patches of the sky. From left to right: reconstructed using Planck-like channels, recon- 
structed using ACT-like channels, reconstructed using jointly Planck-like/ ACT-like channels. From top to bottom: sample of 30 extended 
clusters with 15 arcmin < #200 < 40 arcmin, sample of 30 "typical" clusters with #200 ~ 5 arcmin, sample of 30 compact clusters with 
#200 ~ 3 arcmin. 



MSE (extended) 


Planck- like 


ACT-like 


Joint 




3.95e-ll 


1.23e-10 


3.35e-ll 


MSE ("typical") 


Planck- like 


ACT-like 


Joint 




1.27e-9 


3.74e-10 


8.12e-ll 


MSE (compact) 


Planck- like 


ACT-like 


Joint 




6.52e-7 


3.65e-10 


1.45e-10 



Table 4. Performance of the mean sample profile reconstruction 
measured by the criterion Eq. ( |13| ). 



the ACT-like reconstruction slightly more contaminated due 
to the lack of channels used to clean large scale Galactic fore- 
grounds. Conversely, results in Tab. [4] confirm that Planck- 
like NILC completely fails in reconstructing the mean TSZ 
profile from more compact clusters (#200 < 5 arcmin), due to 
the limited resolution of the channels compared to ACT-like 
channels, whereas the joint Planck-like/ ACT-like NILC sig- 
nificantly improves the profile reconstruction by two orders 
of magnitude. The best bias- variance trade-off is still accom- 
plished by the joint Planck- like/ ACT- like reconstruction of 
the mean TSZ profile. 



where Y t = ±£? =1 1?, = I £» =1 y. 

&i — ~2 5_}j=i {^l ) • Here n is the size of the sample and the 
quantities indexed by j refer to one cluster of the sample. In 
Tab. [4]we have listed the results obtained for the three sets of 
cluster samples (extended, "typical", and compact) and for 
the three sets of exploited data (Planck-like, ACT-like, and 
joint Planck- like/ ACT-like). The results show that the mean 
TSZ profile from extended clusters (#200 > 15 arcmin) is re- 
constructed with similar accuracy for the three datasets with 



and 



4 DISCUSSION 

4.1 Impact of the size of patch and of the number 
of channels 

By contrast with a full-sky ILC analysis, when performing 
component separation on small patches of the sky we en- 
counter the problem of ILC bias discussed in Delabrouille 
et al. (2009). On small patches the statistics is computed on 
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Figure 11. Same as Fig. |10| in logarithmic scale. 
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Figure 12. Residual dust mean profile (extended cluster sample). Planck-like (dot-dashed black), ACT-like (long dashed red), Planck- 
like/ ACT-like combined (solid blue). 
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Figure 13. Impact on NILC reconstruction of the size of the patch and of the number of exploited channels. The reconstruction is shown 
on the selected cluster 1 of Tab.[2]and is performed at a resolution of 10 arcmin. Left panels: the patch is 7° x 7° and all Planck-like and 
ACT-like channels have been exploited by NILC. Middle panels: the patch is 3° x 3° and all Planck-like and ACT-like channels have 
been exploited by NILC. Here a typical ILC bias is observed due to the decreasing of the number of modes in such a small patch of the 
sky. Right panels: in the case of a 3° x 3° patch the bias is removed by decreasing the number of exploited channels to six channels: 3 
Planck-like channels (30 GHz, 90 GHz, and 350) and 3 ACT-like channels (148 GHz, 219 GHz, and 277 GHz). 



a reduced number iV p of mode^jso that it creates artificial 
correlations where the ensemble average would vanish. The 
artificial correlations manifest themselves as a bias in the 
variance of the ILC estimate s~: 



(14) 



where the third term in the right-hand side of Eq. ( 14 ) would 
vanish in the absence of bias, given the assumption that the 
sky components are physically independent. At first order 
in 1/iVp, this bias is an anticorrelation between the signal 
and the residual noise and is given by the simple formula: 



(8(8-8)) 



N D 



(15) 



where iV c h is the number of channels and N p is the num- 
ber of modes used to compute the coefficients of the local 

3 The number of independent modes corresponds to the number 
of needlet coefficients at a given scale (set of neighbouring pixels) 
used to compute the local covariance matrix Eq. |8l. 



covariance matrix Eq. (|8]). This unphysical anticorrelation 
may thus induce a power loss in the TSZ signal reconstruc- 
tion. The exact derivation of formula Eq. (p~5j) can be found 
in the appendix of Delabrouille et al. ( 2009| . It can be in- 
tuitively understood as follows. On the one hand the vari- 
ance of each mode of the true signal is {s 2 )/N p since the 
sum over the modes must have total variance (s 2 ). On the 
other hand, the number of modes available for reconstruct- 
ing the TSZ signal is (N p - (N ch - 1)) because (iV c h - 1) 
degrees of freedom are required to adjust the iV c h weights 
Wi of the constrained minimization problem Eq. ([3| . There- 
fore, the variance of each mode of the reconstructed signal 
by ILC is (s 2 )/(iVp - (iV ch - 1)). The correlation of the re- 
constructed TSZ map with the original TSZ map thus is 
(iV p - (iV c h - 1))/N P . It is this loss of iV c h - 1 modes of the 
original TSZ which induces the negative bias Eq. (15). 



On infinite samples or large patches, iV p goes to in- 



finity so that the bias Eq. (15) goes to zero. Conversely, 



decreasing the size of the patches results in the increase of 



the bias. This is shown in Fig. 13 (middle versus left pan 
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els) where we compare the NILC TSZ reconstruction on a 
3° x 3° patch with the reconstruction on a 7° x 7° patch 
centred at the same location on the sky. In order to reduce 
the bias on small patches we should reduce the number of 
channels used in the ILC filtering, according to formula Eq. 
(15). However, the number of channels has to be significant 
enough in order to enable ILC to correctly filter the con- 
tamination from foregrounds. Therefore, we are left with a 
compromise between minimizing the signal bias versus min- 
imizing the background noise. On the right panels of Fig. 
J3]we show that we are able to remove the ILC bias on a 
3° x 3° patch TSZ reconstruction by reducing the number 
of exploited channels to six (30 GHz, 90 GHz, and 350 GHz 
Planck-like channels, and the three ACT-like channels) while 
this number of channels is kept significant to minimize the 
background noise from the other sky components. We delib- 
erately selected a subset of six channels covering the whole 
range of available frequencies going from 30 GHz to 350 GHz 
in order to guarantee a robust cleaning of the residuals from 
other sky emissions in the NILC TSZ reconstruction. 



4.2 Type of wavelets 

We have shown in Sect. 12.31 that needlets enables ILC to 
combine heterogeneous instrument datasets with different 
resolutions for TSZ reconstruction. For comparison, we re- 
peat TSZ component separation combining both Planck-like 
and ACT-like observations with an ILC based on different 
types of wavelets: Daubechies wavelets, symlets, coiflets (see 
Daubechies (1992) and the references therein) and needlets 
(Guilloux, Fay Sz Cardoso[|2009|. The definition of the stan- 



dard wavelets (Daubechies, symlets, coiflets) is such that the 
dilation factor is sampled on a dyadic scale. The results on 
the reconstruction of the TSZ signal from the selected clus- 
ter 1 of Tab. [2] are shown in Figs. [l4| and |15[ 

We display in Fig. [14] the scattering of the wavelet-ILC 
reconstructed signal with respect to the input TSZ signal. 
The wavelets showing the minimum residual noise are the 
symlets of order 8, the coiflets of order 5, and the needlets. 
Among these three types of wavelet-ILC, the needlet ILC 
shows a better reconstruction of the cluster profile in Fig. 
\TE\ the other two being slightly biased. The robustness of 
the wavelet-ILC reconstruction depends on the properties 
of localization of each type of wavelet both in space and 
in scale. The fixed dyadic sampling of the Fourier scales 
probed by the standard wavelet-ILC is responsible for a re- 
duced number of modes in the reconstruction of the signal 
whereas a needlet ILC appears more flexible to sample the 
Fourier domain (number and width of needlet bands freely 
tunable) as discussed in |Guilloux, Fay &; Cardoso] (12009 ) . 
The loss of Fourier modes by a fixed sampling of the wavelets 
may induce some bias in the reconstruction compared to the 
needlet ILC reconstruction. From the basic tests made here, 
we can see in Figs. [l4| and [T5] that the best compromise be- 
tween a minimum bias on the TSZ profile and a minimum 
background noise in the reconstructed signal is accomplished 
by the use of needlets. A detailed analysis on wavelet local- 
ization would be required to optimize the TSZ component 
separation. 



5 CONCLUSION 

The reconstruction of galaxy cluster profiles through ther- 
mal SZ effect is an exciting challenge in data analysis of cur- 
rent CMB experiments because the SZ profile is a powerful 
probe of the baryon physics in the cluster outskirts, beyond 
the radius of virialisation in the intra-cluster medium. In this 
region, X-ray measurements fail to reconstruct the thermal 
pressure profile. For the first time, the SZ profiles of some 
galaxy clusters have been reconstructed from the surveys of 
the high resolution telescopes ACT ( |Sehgal et al. 2011 ) and 
SPT ( |Plagge et ah] |2010| ). Results on the SZ profile from 
the Planck survey have also been published very recently 
in |Planck Collaboration et al.| ( |2012a| ). Both in ACT and 
SPT results, the SZ signal has been reconstructed by us- 
ing a matched spatial filtering, which is an effective way of 
reconstructing the SZ effect when few frequency maps are 
available from the observations. However, matched filtering 
relies on priors on the expected template profile, which may 
appear controversial when one is interested in reconstructing 
the cluster profile. The needlet ILC on patches developed in 
this work is a blind approach which does not suffer from any 
prior to reconstruct the SZ profiles. 

Apart from the prior issue, the accuracy of the SZ re- 
construction also relies on the achievable resolution of the 
instrument and on the level of contamination either by sky 
emissions (Galactic foregrounds, CMB, etc) or by the in- 
strumental noise. We have highlighted that both problems 
can be tackled simultaneously by combining extra frequency 
maps from multiple instrument datasets with different res- 
olutions. We have shown in this work how needlet ILC has 
the ability to adapt the needlet windows in Fourier space 
to the beams of the channel-maps exploited for component 
separation, thus offering the possibility of combining het- 
erogeneous datasets coming from multiple instruments with 
various resolutions and different frequency coverages. On 
the one hand, the properties of localization of the needlets 
both in scale and in space enable component separation to 
adapt to the local conditions of contamination (physical at 
large scales, instrumental at small scales), on the other hand 
they enable component separation to aggregate heteroge- 
neous instrument datasets with different achievable resolu- 
tions. We have applied needlet ILC on three different simu- 
lated datasets (single Planck-like, single ACT-like, and joint 
Planck- like/ ACT- like) for different patches of the sky cen- 
tered on various extended and compact clusters. The per- 
formance of the reconstruction has been validated both on 
particular clusters and on samples of selected clusters. The 
aggregation of multiple datasets with needlet ILC allows us 
to reconstruct the SZ profile over a large radius with both 
high resolution and robust foreground and noise cleaning. 

We have also explored the limitations of the method. 
Needlet ILC applied on a small-size patch of the sky encoun- 
ters a bias in the reconstructed SZ signal. However, combin- 
ing a subset of channels in that case enables ILC to tackle 
this problem of bias. An optimized selection of the subset 
of channels would be required to improve the reconstruction 
on very small patches. It would also be very instructive to 
optimize in a future work the wavelet localization (type of 
wavelets, width of spectral windows) for multi-instrument 
SZ component separation. 

Needlet ILC performed on the joint dataset in- 
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Figure 14. Wavelet ILC comparison (scatter plots) for the selected cluster 1 of Tab. [2] TOP : Daubechies wavelets of order 2, 8, and 24. 
MIDDLE : Symlets of order 2, 4, and 8. BOTTOM : Coiflets of order 1, 5, and Needlets. 



eluding both Planck-like and ACT-like datasets success- 
fully reconstructs the TSZ effect from compact clusters 
(#200 < 5 arcmin) beyond the beam of the single Planck-like 
instrument while cleaning foreground residuals better than 
a component separation applied on the single ACT-like in- 
strument. The multiresolution approach presented in this 
work appears to be promising for a future multi-instrument 
analysis of the SZ effect. 
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Figure 15. Wavelet ILC comparison (TSZ profile) for the selected cluster 1 of Tab. [2] TOP : Daubechies wavelets of order 2, 8, and 24. 
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